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Emergent spin-1 trimerized valence bond crystal in the spin-1/2 Heisenberg model on the star lattice 
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We explore the frustrated spin-1/2 Heisenberg model on the star lattice with antiferromagnetic (AF) coup¬ 
lings inside each triangle and ferromagnetic (FM) inter-triangle couplings (Jg < 0), and calculate its magnetic 
and thermodynamic properties. We show that the FM couplings do not sabotage the magnetic disordering of 
the ground state due to the frustration from the AF interactions inside each triangle, but trigger a fully gapped 
inversion-symmetry-breaking trimerized valence bond crystal (TVBC) with emergent spin-1 degrees of free¬ 
dom. We discover that with strengthening Jg, the system scales exponentially, either with or without a magnetic 
field h: the order parameter, the five critical fields that separate the Je-h ground-state phase diagram into six 
phases, and the excitation gap obtained by low-temperature specific heat, all depend exponentially on Jg. We 
calculate the temperature dependence of the specific heat, which can be directly compared with future experi¬ 
ments. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.60.Ej, 05.70.Ln 

Introduction .— Two-dimensional (2D) spin-1/2 frustrated 
magnetic systems [1] are currently of great interest, because 
they may realize exotic quantum states that do not possess any 
semi-classical spin ordering [2], such as quantum spin liquids 
(QSLs) or valence bond crystals (VBCs). Leading candid¬ 
ates for realizing such states are spin-^ Heisenberg models 
with competing interactions on, e.g. square, honeycomb and 
kagome lattices [3-17]. A particularly promising QSL system 
that has been argued to have experimental realizations is the 
kagome Heisenberg antiferromagnet (KHAF) [18-24]. How¬ 
ever the nature of its ground state, i.e., a gapped Z 2 spin liquid 
[3-6] versus a gapless U(l) Dirac spin liquid [7-10], is still 
under debate. 

Another frustrated 2D quantum system of great potential 
interest is the Heisenberg model on a star lattice (Fig. 1). Its 
physics is arguably even richer than that of the KHAF, for sev¬ 
eral reasons: (a) similar to the kagome lattice, the star lattice 
bears a high geometrical frustration due to its triangle struc¬ 
ture; (b) the star lattice possesses a lower coordination num¬ 
ber than the kagome lattice, implying stronger fluctuations; 

(c) the star lattice naturally involves two inequivalent bonds, 
which can lead to exotic quantum phases; (d) various QSLs, 
such as the non-Abelian chiral spin liquid and the double se- 
mion spin liquid, have been found in several models on a star 
lattice, e.g., the Kitaev model and the quantum dimer model 
[25, 26]; (e) a number of organic iron acetates have been syn¬ 
thesized in experiments [27], which can be described by the 
Heisenberg model on a star lattice. 

However, the properties of the Heisenberg model on a star 
lattice have not been fully explored yet. Recent research using 
the large-W approximation and a Gutzwiller projected wave- 
function [28] only investigated the system for antiferromag¬ 
netic (AF) inter-triangle couplings (Je > 0), where a Je-dimer 


VBC and a \/3 x \/3 VBC phase [29] are found (Fig. 1). How¬ 
ever, the properties of the system for the ferromagnetic (FM) 
Je < 0 are still unclear, where the system can be closely re¬ 
lated to spin-1 kagome physics [30-32]. 

The intrinsic importance of 2D frustrated many-body sys¬ 
tems is matched by the great technical challenges involved 
in studying them. One such challenge is calculating thermo¬ 
dynamic properties, such as speciflc heat and susceptibility, 
which can be detected directly in experiments. Most existing 
studies have, however, focused on the ground-state. Indeed, 
numerical studies of frustrated quantum magnetic systems at 
flnite temperature are very scarce, owing to the extreme diffi¬ 
culties of such calculations [33-35]. 

In this Letter, we employ four different state-of-the-art al¬ 
gorithms [34, 36-41] to perform a comprehensive study of the 
spin-^ Heisenberg antiferromagnet on the star lattice with FM 
inter-triangle couplings (Je < 0), calculating its ground state 
and thermodynamic properties. We show that the FM inter¬ 
triangle couplings do not sabotage the magnetic disordering of 
the ground state that arises due to frustration generated by AF 
intra-triangle couplings, but, remarkably, trigger a trimerized 
valence bond crystal (TVBC) with emergent spin-1 degrees of 
freedom, that breaks spatial inversion symmetry. We determ¬ 
ine the phase diagram of the system in a magnetic field and 
identify six phases. We uncover a magnetization cusp on the 
boundary between the inversion-symmetry-breaking and the 
non-inversion-symmetry-breaking phases. We calculate the 
temperature dependence of the speciflc heat and determine a 
non-magnetic gap by analyzing accurate results for the low- 
temperature behavior of the speciflc heat. We And that the 
observables of the system scale exponentially with increas¬ 
ing I Je|, evidenced by the large-Je dependence of a range of 
physical quantities, such as the TVBC “order parameter”, five 
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Figure 1: (Color online) The ground-state phase diagram of the star 
Heisenberg model. For Je > 0, previous studies show various pos¬ 
sible VBCs and spin liquids, where one recent work found a \/3 x \/3 
VBC and a Je-bond VBC [28, 29]. The phase boundary has not 
been settled yet. For Je < 0, we show that the system is in a trimer- 
ized valence bond crystal (TVBC) phase, where a triplet appears at 
each Je bond and the inversion symmetry of up and down triangles 
(marked by blue and yellow, respectively) is broken. 

critical fields and the non-magnetic gap. 

Model and methods .— The Hamiltonian of the star Heisen¬ 
berg model reads 



Figure 2: (Color online) (a) The ground-state energy Eq of the star 
Heisenberg model for Je = — 1 and Jt = 1 obtained by the plain and 
SU(2) simple update algorithm, where Eq converges versus number 
of states {D) and number of multiplets (D* red solid dots), respect¬ 
ively, with clearly superior performance of the SU(2)-based calcu¬ 
lations. For comparison only, we also translate the number of mul¬ 
tiplets D* into the corresponding actual number of states D (black 
lines and symbols). Since for fixed D*, D may vary (depending 
on the specific set of multiplets associated with each virtual bond), 
we show the minimum, maximum and average value of D over the 
three virtual bond indices of a tensor, (b) Eq versus Je, obtained by 
SU(2) simple update and DMRG simulations, which show very good 
agreement with each other in the whole parameter range. The inset 
sketches the local tensors of the TN state. 


H = Je E Si-Sj + Jt E 

{ij)eJe {im)eJt 

where the first summation runs over all inter-triangle bonds 
and the second over all intra-triangle bonds. We employ 
tensor network [37] and DMRG [39] methods to simulate the 
model on the infinite lattice and cylindrical geometries, re¬ 
spectively. To be specific, the TN representation of the ground 
state [inset of fig. 2(a)] can be written as 


|V’) = ^Tr{,}eTN[n™ 

{«} i 


Sj,lSj,2Sj,3 




where T{j) is a (d^ x D^) tensor residing on the j-th tri¬ 
angle with physical dimension d and ancillary bond dimen¬ 
sion D, containing all parameters of the TN state. The an¬ 
cillary bonds {aj^n} (e = 2, 3) carry the entanglement of 

the state and Tr|a}eTN denotes a contraction of all shared 
{cij^n}‘ The physical bonds = 1, 2, 3) represent 

the three spins inside the j-th triangle with local basis \sj^n)^ 
Such a TN ansatz is called a projected entangled-pair state 
(PEPS) [40]. The simple update algorithm provides an effi¬ 
cient way to optimize the PEPS by minimizing the energy per 
siteJ;o = (t/;|JJ|t/;). 

In addition, we implement SU(2) symmetry in TN states 
and related algorithms by using QSpace techniques [41]: 
we impose SU(2) symmetry in every single tensor index, 
retain the symmetry during imaginary time evolutions and 
other tensor manipulations, keep track of multiplets (instead 
of individual states) on the bonds, and manipulate only re¬ 
duced tensors (instead of full tensors), thus reducing both the 
memory and CPU time dramatically [30, 41]. 


Eirst, we compare the ground-state energy obtained by dif¬ 
ferent methods. In Eig. 2(a), we show the energy obtained by 
plain and SU(2) PEPS calculations, which both converge to 
the same results. Note that for comparable number of states 
J9, a lower ground state energy can be obtained by plain PEPS 
as compared to SU(2) PEPS, since it is allowed to break sym¬ 
metries and hence has access to a larger variational parameter 
space. However, the results converge towards the same value 
for large D, suggeting that as expected, the tensors eventually 
converge to tensors that respect symmetries. This justifies the 
exploitation of symmetries at significantly reduced overall nu¬ 
merical cost. 

In Pig. 2(b) we show the energy obtained from SU(2) TN 
simulations and cylindrical DMRG for —10.0 < Je < 3.0, 
which show an excellent agreement in the whole region. 
Moreover, the appearance of a cusp in the energy curve at 
Je = 0 indicates a first-order phase transition. 

Spontaneous inversion symmetry breaking .— We now 
study the ground state of the star model, which is found to 
possess spontaneous inversion symmetry breaking (SISB). It 
can be characterized by the energy difference between the two 
kinds of triangles S = \E^ — E^\ where we have = 

{^\ per triangle with the summation run¬ 

ning over all local interactions Hij inside the up (down) tri¬ 
angles. We use DMRG to calculate the cylinder system with 
the geometry shown in Pig. 3(a) (denoted by YC4). To break 
the inversion symmetry between the up and down triangles, 
we take the couplings inside the up triangles on the open 
boundaries as Jpin = 2Jt (Jt is the coupling constant for all 
other triangles). Then, we can measure the decaying behavior 
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(a) Bond correlation, Je=-4, YC4-36 cylinder 



Figure 3: (Color online) (a) The bond correlation (Si • Sj) versus 
distance x from the pinned boundary (using Jpin = 2) on a YC4- 
36 cylinder for Je = — 4 calculated by DMRG keeping 2000 SU(2) 
multiplets. (b) Log-linear plot of the inversion symmetry breaking 
parameter S as function of the distance x from the boundary for the 
YC4-36 cylinder with boundary pinning, (c) The Je dependence of S 
obtained from SU(2) simple update simulations. As long as Je < 0, 
the system has a non-zero S. We find by fitting that for approximately 
— Je ^ 4, (5 fulfills the relation S = S{1 — where = 0.28 

and (5 = 0.1 that gives the value of (5 for Je ^ — oo. In contrast, for 
small I Je| we find that S satisfies S = 0.03Je for Je ^ 0 (see inset). 

of S from the boundary to the bulk. As shown in Fig. 3(b), we 
find that S decays quite slowly, implying a large decay length 
for large | Jg |. We checked that different values of Jpin give 
the same decay length. Since a non-symmetry-breaking state 
just possesses a short decay length, our calculations strongly 
suggest a ground state with SISB on a wider system. With 
decreasing | Jg |, S decays faster. For small values of | Jg |, the 
SISB is too weak to identify on a small cylinder. Thus, we 
employ the TN simulations on the infinite-length system. 

The order parameter S obtained from TN methods is shown 
in Fig. 3(c). For large | Jg|, TN calculations, too, find a strong 
TVBC order, consistent with DMRG results. By fitting S with 
—Jg 0, we find that S fulfills an exponential behavior with 
Jg as 

S = S(1 - (3) 

where we have fi = 0.28 and 5 = 0.1. It indicates that the 
large | Jg| couplings project each corresponding spin-^ pair 
into an effective S = 1 spin, and stabilize a TVBC. Inter¬ 
estingly, for the small | Jg| region, the TN simulations show 
that the inversion symmetry is broken for any small Jg < 0, 
while such a symmetry is found to be intact for Jg > 0. To 


be specific, for — Jg —> 0, (5 satisfies the algebraic relation 
5 = 0.03Jg, as shown in the inset of Fig. 3 (c). Our results 
not only support the TVBC ground state for the spin-1 kagome 
model [30-32], but also further show that such a TVBC is ro¬ 
bust in the spin-^ star model for any finite FM Jg interaction, 
where each two spin-^ connected by Jg cannot be regarded as 
strictly projected into an effective spin-1. 

Ground-state phase diagram in magnetic fields. — In a 
magnetic field, frustrated magnetic systems usually exhibit 
singularities in the magnetization curve such as cusps [42] and 
plateaus [43], which reveal the exotic structure of the energy 
spectrum and distinguish different phases. We study the field 
dependence of the magnetization per site Mz and the energy 
difference J, as shown in Fig. 4. Interestingly, we find a zero 
plateau corresponding to a finite spin gap, a cusp represent¬ 
ing the restoration of inversion symmetry, and a 1 /3-plateau 
in the magnetization curve. 

In the zero plateau region, h < hd , both Mz and 5 remain 
unchanged, indicating that there is a finite spin gap protecting 
the TVBC state. With increasing h, the spin gap decreases 
and eventually closes at = hd . For h > hd , becomes 
nonzero and 5 starts to decrease. Ath = hc 2 , a cusp appears in 
Mz and 6 vanishes, separating the SISB phase from non-SISB 
phases. A magnetization cusp has also been observed in some 
one-dimensional frustrated magnetic systems having ground 
states that break lattice symmetry, reflecting the novel energy 
dispersion of the low-lying excitations [42]. A first shoulder 
in the magnetization occurs consistently around M 1/30. 
By further increasing the field, we And a 1 /3-plateau corres¬ 
ponding to a gapped solid state [43]. Based on the behaviors 
of Mz and 6 we obtain the quantum phase diagram in the Je-h 
plane, shown in Fig. 4(c). 

We And that the critical fields hd (i = 1, 2, • • • , 5) also 
converge exponentially for large | Je |, 

hd = hd{l - (4) 

as shown in Fig. 4(c), with coefficients given in Table I. The 
scaling behavior of the critical fields strongly implies that the 
star Heisenberg model approaches the effective spin-1 model 
in an exponential manner, suggesting that the large | Jg | rep¬ 
resents a gapped system, consistent with [30-32]. 

Specific heat. —To study the thermodynamics and the low- 
lying excitations, we calculate the temperature dependence of 
the specific heat, employing the network contractor dynamics 
(NCD) method [34] that is based on a TN representation of the 
finite-temperature density matrix. In experiments, the specific 
heat can be measured by mature techniques, e.g. a thermal 
relaxation calorimeter. 

As shown in Fig. 5(a), by taking Jg from zero to — oo, 
the low-temperature peak moves to higher temperature and 
merges with other peaks. From the inset of Fig. 5(a), one can 
see that below the low-temperature peak. In C depends lin¬ 
early on the inverse temperature 1/T aslnC = —AjT + 
const., indicating a finite gap A that is consistent with the 
gapped TVBC ground state. The Jg-dependence of A is given 
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Figure 4: (Color online) The field-dependence of (a) the magnetiza¬ 
tion Mz and (b) the TVBC order parameter 5. Five critical fields hd 
{i — 1, 2, • • • , 5) is determined by Mz and S, which determine six 
phases in the Je-h diagram as shown in (c). For 0 < h < hd, 
Mz — 0 and 5 is intact. For hd < h < hc2, Mz increases, 
and S starts to diminish and vanishes at = hc 2 , where the inver¬ 
sion symmetric and broke phases are separated and one always has 
Mz = 1/30. For hc3 < h < hc4, the system is in a conventional 
1 /3-plateau solid phase. 

in Fig. 5(b). We observe that A also converges exponentially 
for large | Je| 

A = (5) 

where a fit yields k = O.b and A = 0.17±0.02 correspond¬ 
ing to the gap for Je ^ — cxo. We suggest that the A obtained 
via the specific heat is the non-magnetic singlet-singlet (not 
singlet-triplet) gap, since SU(2) symmetry is preserved at all 
temperatures, evidenced by the vanishing finite-temperature 
magnetization. This is consistent with the spin-1 kagome cal¬ 
culation [44], where a non-magnetic gap A = 0.1 ^ 0.2 is 
obtained. 

The non-magnetic gap is larger than the spin gap, imply¬ 
ing that there are no singlet states lying under the first triplet 
excitation, thus the TVBC is fully gapped. This is consistent 
with the fact that the magnetic field hc 2 , at which inversion 
symmetry is recovered, is larger than hd, at which the spin 
gap closes. 

Conclusions. — We discover an emergent spin-1 TVBC 
with spontaneous lattice-inversion-symmetry breaking in the 
spin-^ star Heisenberg model with FM inter-triangle coupling 
Je, and study its ground-state and thermodynamic properties. 
The reliability is confirmed by four different algorithms in¬ 
cluding SU(2) DMRG, simple update of the TN state with and 
without SU(2) symmetry, and NCD. Our results not only sup¬ 
port the conclusion of a TVBC state in the ground state of the 
spin-1 kagome Heisenberg model [30-32, 44], but also reveal 
the rich properties of the spin-^ star system, including fruit¬ 


Table I: Values for the fitting parameters hd, m and Ui of the critical 
fields [see Eq. (4)]. Note that as hd = 1.5 is a constant, we have 
05 = 0 and any 



hci 

hc2 

hc3 

hcA 

hc5 

hci 

0.07 

0.132 

0.34 

0.81 

1.5 

OLi 

1 

1 

1 

-0.85 

0 

Ui 

0.5 

0.65 

0.7 

0.6 

* 


T 

Figure 5: (Color online) (a) The temperature-dependence of specific 
heat C for various Je, where multi-peak structures are observed. The 
position of the low-temperature peak moves to the higher temperat¬ 
ure as — Je increases. Inset: the curve of In C versus inverse temper¬ 
ature 1/T. Below the low-temperature peak, one can see that the spe¬ 
cific heat decays exponentially with 1/T as In C = — A/T + const, 
where A is the (Je-dependent) excitation gap. (b) By fitting specific 
heat, the excitation gap A for different Je is obtained and shown to 
fulfill the relation A = A(1 — e'^'^®). The error bars are given by the 
linearity of the low temperature C. By fitting the A — Je curve, we 
have A = 0.175 that gives the excitation gap in the Je ^ — oo limit 
and the constant n = 0.5. 


ful phases in a magnetic field, magnetic cusps at Mz ^1/30 
and exponential scaling behavior for Je — oc. Our work 
implies that spin-1 VBCs may emerge in the geometrically 
frustrated spin-^ systems with FM interactions. Moreover, 
our calculations of the specific heat present an efficient way 
to obtain the excitation gap of 2D frustrated many-body sys¬ 
tems. It also provides valuable information at finite temperat¬ 
ures which can be compared directly with future experiments. 
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